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ABSTRACT 
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The  height  dependence  of  the  temperature  structure  parameter,  C  ' 
has  been  measured  with  microthermal  sensors  mounted  on  a  light  air- 
craft. This  work  was  done  in  conjunction  with  optical  propagation 
and  turbulent  transport  research  in  the  marine  boundary  layer. 
These  measurements  indicate  that,  in  the  absence  of  a  strong  inver- 
sion, the  constant  stress  layer  can  be  surprisingly  thin.  The  measure- 
ments also  substantiate  the  strong  role  played  by  temperature  and  water 
vapor  discontinuities  in  turbulence  above  the  boundary  layer. 
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I .   INTRODUCTION 

2 
We  have  made  measurements  of  temperature  structure  parameter,  CL  ,  from 

a  light  aircraft  using  microthermal  sensors  as  part  of  a  study  of  turbulence 

2 
in  the  marine  boundary  layer.  CT  is  important  for  optical  propagation 

studies  due  to  its  relation  to  the  index  of  refraction  structure  parameter, 

CN2,  [Friehe  (1977)] 

CN2  =  (79  x  10"6P/T2)2  (CT2  +  .11  CL  +  3.2  x  10"3C  2)  (1) 

2 
where:    C    is  the  water  vapor  structure  parameter, 

CT   is  the  cospectrum  structure  parameter, 

p    is  the  pressure  in  mb,  and 

T    is  the  absolute  temperature. 

The  water  vapor  fluctuations  are  usually  relatively  small  so  we  can  write 

C^2  =  (79  x  10"6P/T2)2  CT2  (2) 

2   2 
This  relationship  is  shown  in  Fig.  1  with  C^  /CL,  as  a  function  of  altitude 

for  the  U.S.  Standard  Atmosphere.  Turbulence  is  also  of  interest  because  of 
the  role  of  eddy  diffusion  in  the  transport  of  heat,  water  vapor,  and  pol- 
lutants. These  factors  are  important  in  the  formation  of  marine  fog  and  air 
pollution  modeling. 

2 
Only  a  few  measurements  of  the  altitude  dependence  of  CL,  have  been  made 

2 
to  date.   Korpov  and  Tsvang  (1966)  have  made  measurements  of  CL  with  an 

acoustic  anemometer  on  an  airplane  and  have  related  their  results  to  the 

vertical  temperature  gradient.  Microthermal  sensors  have  been  used  to  measure 

2 
CL  with  balloon  borne  equipment  by  Bufton  (1973)  and  airplane  borne  equipment 


by  Lawrence  (1970)  and  Collins  (1977) .  Recently,  Hanson  (1976)  has  combined 

airplane  microthermal  measurements  and  remote  scintillometer  determinations 

2 
of  C,  and  compared  his  results  with  an  impirical  model  developed  by  Yura. 

Based  on  this  body  of  data,  Hall  (1977)  has  compiled  a  conglomerate  curve 

2 
of  C_  as  a  function  of  altitude  (Fig.  2)  for  daytime  overland  data.  Only 

Ochs  (1973)  has  reported  measurements  over  the  ocean. 

The  data  we  are  reporting  on  were  taken  as  part  of  three  separate 
research  operations .  The  first  was  the  "Cooperative  Experiment  for  West 
Coast  Oceanography  and  Meteorology  -  1976"  (designated  CEWCOM  76) ,  a  marine 
fog  research  project  off  San  Diego.   CEWCOM  76  was  organized  by  the  Naval 
Postgraduate  School  and  the  Naval  Electronics  Laboratory  Center.  The  second 
was  a  marine  fog  and  aerosol  project  in  the  Gulf  of  Mexico  off  Panama  City, 
Florida  (designated  FLORIDA  77) .  The  third  was  in  conjunction  with  a  turbu- 
lence and  aerosol  research  cruise  on  the  USNS  Hayes  in  the  Atlantic  (desig- 
nated HAYES  77)  organized  by  the  Naval  Research  Laboratory. 


II.   INSTRUMENTATION  AND  TECHNIQUES 

The  platform  for  our  measurements  is  a  single  engine  turbocharged 
Bellanca  operated  by  Airborne  Research  Associates  (Fig.  3) .  This  aircraft 
has  been  flown  as  low  as  3  meters  and  as  high  as  10,000  meters  and  makes  an 
excellent  tool  for  the  low  altitude  flights  required  for  boundary  layer 
research.  The  aircraft  is  well  instrumented,  allowing  simultaneous  measure- 
ments of  air  temperature,  altitude,  dew  point,  electric  field,  visibility, 
infrared  surface  temperature,  and  microwave  refractive  index.  The  data  is 
normally  recorded  with  an  eight  channel  strip  chart  recorder. 

The  temperature  structure  parameter  is  measured  using  the  paired  sensor 
method.   Given  two  temperature  sensors,  1  and  2,  a  distance  d  apart,  then 
[Lumley  (1964)] 

CT2  =  <  (T2  -  T:)2  >  d'2/3  (3) 

This  quantity  is  related  to  the  Kolmogorov  power  spectral  density  of  tempera- 
ture fluctuations,  <j>T(k), 

*TCk)  =  .25  CT2  k~5/3  (4) 


where  k  is  the  wave  number.  Equation  4  applies  in  the  inertial  subrange  part 

of  the  spectrum  where  the  turbulence  is  nearly  isotropic,  allowing  a  one- 

2 
dimensional  representation  of  <J>T(k)  .   C™  is  independent  of  d  in  the  inertial 

subrange . 

2 
The  device  we  have  used  to  measure  CL,  is  a  DC  Wheatstone  bridge  (Thermo 

Systems  Model  1044)  that  senses  the  relative  resistance  fluctuations  of  a 


pair  of  2.5  micron  diameter  platinum  wires  separated  by  about  one  meter. 
The  sensors  were  originally  mounted  on  tne  leading  edge  of  the  wing  (CEW- 
COM  76  and  FLORIDA  77)  but  the  noise  level  was  very  high  so  the  sensors 
were  moved  to  the  present  wing  tip  location  (Fig.  4) .  The  output  of  the 
bridge  (proportional  to  AT  =  T2  -  T,)  is  processed  by  an  RMS  module  with 
a  5  sec  time  constant  and  recorded  on  the  strip  chart  recorder  as  ATRMC. . 
The  sensitivity  is  limited  by  the  broad  band  noise  of  the  system  due  to 

inherent  amplifier  noise  and  pickup  from  the  aircraft  ignition  system. 

2        -4 
For  CEWCOM  76  and  FLORIDA  77  the  noise  level  corresponds  to  L  ~  4  x  10 

2  2/3  2     -4  2  2/3 

k  /m    ,  for  HAYES  77  the  noise  level  is  lower  (CL  — '  10   k  /m   )  since 

we  were  using  the  wing  tip  probe  configuration.   It  is  possible  to  improve 

the  accuracy  by  correcting  for  the  noise.   If  we  assume  an  RMS  noise  level 

of  N,  then  the  noise  correction  appears  as 

ct2  =  c  <AW2  -  n2  i  d"2/3  & 

Since  2.5  micron  wires  are  fragile,  breakage  is  a  continual  problem. 
Except  for  one  bad  batch  of  wires,  we  have  found  a  typical  lifetime  of  one 
or  two  flights  for  a  given  wire.  We  presently  have  two  pairs  of  sensors 
mounted  and  can  switch  to  a  good  pair  if  one  wire  breaks . 


III.   THEORY 

The  boundary  layer  is  that  part  of  the  atmosphere  where  friction 
with  and  heating  by  the  surface  play  an  important  part  in  the  generation 
of  turbulence.  Near  the  surface  the  shear  stress  and  scalar  fluxes  are 
essentially  constant.   In  this  region  the  fluxes  can  be  represented  by 
scaling  parameters  (such  as  U#  and  T*)  that  are  independent  of  height. 
We  shall  refer  to  this  layer  of  nearly  constant  stress  and  flux  as  the 
surface  layer.   In  the  surface  layer,  the  height  above  the  surface,  Z, 
is  the  appropriate  turbulence  length  parameter.  For  a  complete  treatment 
of  the  surface  layer  equations  we  suggest  Lumely  (1964)  ,  Businger  (1971) 
or  Kraus  (1972) .  The  normalized  momentum  flux,  F  ,  and  the  normalized 
heat  flux,  F,  ,  for  turbulent  transport  are 

F  =  -  <  u'w1  >  =  LL2  (6) 

m  * 


Fh  =  -  <  T'w'  >  =  U*T*  (7) 


where  u'  is  the  horizontal  velocity  fluctuations, 

w"  is  the  vertical  velocity  fluctuations, 

T'  is  the  temperature  fluctuations, 

U*  is  the  friction  velocity,  and 

T*  is  the  scaling  temperature. 

The  atmospheric  stability  is  represented  by  either  the  Monin-Obukov  length, 

L,  or  the  Richardson  number,  R. , 

l 


2 
T  U* 

L  ■  wz  (8) 


g(3T/3Z) 

R. V—j  (9) 

T(3U/3Z) 


where  T  is  the  virtual  potential  temperature, 

g  is  the  acceleration  of  gravity, 

U  is  the  mean  horizontal  velocity,  and 

K  =  .35  is  the  Von  Karmon  constant. 

The  mean  and  fluctuating  temperature  dependences  on  height  are  given  by 
[Wyngaard  (1971)] 

T 

I  =  kz   fi  w  (10) 


CT2  =  T*2  Z_2/3  f2  (Z/L)  (11) 


Under  near-neutral  conditions  f,  and  f?   are  equal  to  unity,  resulting  in  a 

2  -2/3 

logarithmic  mean  temperature  profile  and  CL  proportional  to  Z  '   .  Under 

2  -4/3 

unstable  conditions  (-  Z/L  »  1/7)  CT  is  proportional  to  Z    .  These 

relationships  are  based  on  measurements  made  on  a  flat  Kansas  plain  with 
averaging  times  of  about  one  hour.  Due  to  the  shorter  average  time  invol- 
ved in  aircraft  profiles,  one  expects  scatter  about  the  curve  of  Eq.  11 
even  in  the  surface  layer. 

Wyngaard  points  out  that  although  these  results  are  based  on  surface 

measurements  (Z  <  22  m)  they  are  valid  to  somewhat  greater  heights.   In 

2 
the  case  of  the  well  developed  unstable  boundary  layer,  the  predicted  CL, 

profile  is  valid  well  beyond  the  surface  layer.  Although  Davidson  (1977) 

has  found  evidence  of  wave  influence  restrictions  on  the  lower  limit  of 


the  oceanic  surface  layer  equations,  the  primary  interest  of  optics  and 
transport  users  is  in  establishing  the  upper  limits  of  validity. 

The  upper  limit  is  often  assumed  to  be  at  least  half  the  distance 
to  the  first  inversion.  This  distinction  becomes  even  more  tenuous  if 
there  is  no  low  inversion.  Above  the  surface  layer,  equations  10  and  11 
become  meaningless  when  defined  in  terms  of  the  absolute  height  above  the 
surface,  Z.   Since  the  vertical  gradients  are  still  meaningful,  Richardson 
number  remains  a  useful  representation  of  stability.  The  appropriate 
length  parameter  is  the  integral  scale  (or  outer  scale),  A,  which  is  a 
measure  of  the  largest  size  (or  minimum  wave  number)  for  which  the 
Kolmogorov  spectrum  of  equation  4  is  valid  [Hinze  (1959)].   In  the  surface 
layer  A  is  proportional  to  Z .   In  this  representation  we  have 

I  =  Fat  fi  (V  ™ 

CT2  =  T*2  (AT)~2/3  ff»  (R.)  (13) 


where  f,'  and  f '  are  analogous  to  f.  and  £~ 


IV .  RESULTS 

2 
Shipboard  and  platform  measurements  of  C_  have  shown  fairly  good 

agreement  with  the  predicted  height  dependence  in  the  near  surface  layer 
(Z  <  25  meters).   Fig.  5  shows  a  profile  taken  at  the  Naval  Coastal  Systems 
Laboratory's  Stage  I  in  the  Gulf  of  Mexico  during  FLORIDA  77.  We  have 
found  the  marine  surface  layer  to  be  predominately  near -neutral  with  a 
tendency  to  be  slightly  unstable.  Based  on  numerous  shipboard  measure- 
ments we  have  found  typical  values  of  L:  •  100  meters  and  T*  r  -  .08  °C. 

2 
Using  these  values  we  have  indicated  typical  surface  based  CT  profile 

(from  equation  11)  as  the  dashed  line  in  Fig.  2.  Under  these  conditions 

we  would  expect  C„  **   Z     to  be  a  good  approximation  for  Z  >  20  meters. 

However,  during  HAYES  77  we  found  the  Atlantic  Coast  from  Cape  Code  to 

Newfoundland  to  have  a  stable  surface  layer.   It  has  been  our  experience 

from  various  shipboard  operations  that  stable  conditions  are  most  likely 

2 
to  produce  anomalous  CT  profiles  in  the  near  surface  layer. 

The  surface  layer  is  usually  well  defined  off  the  Pacific  Coast  of 
the  United  States  due  to  the  persistence  of  a  strong  marine  inversion. 

Consequently,  we  can  expect  CT  to  be  well  described  by  the  Z     or  Z 

2 
equation.   In  Fig.  6  we  have  two  aircraft  measurements  of  CT  before  and 

after  a  radiosonde  balloon  launch  during  near  neutral  conditions.  The 

CT  profile  is  very  well  fit  by  the  Z~    law  until  the  inversion  is 

2 
reached,  where  CT  increases  rapidly  with  the  temperature  gradient.   In 

this  case  the  surface  layer  dominates  the  entire  boundary  layer.  The  strong 

2 
peak  in  CT  at  the  inversion  (Z  ~  200  m)  is  in  agreement  with  the  ship's 

acoustic  sounder. 


The  Atlantic  Coast  data  taken  during  HAYES  77  is  considerably  less 
encouraging.   In  Fig.  7  we  can  see  a  well  developed  surface  layer  similar 


to  the  Pacific  Coast  profile  of  Fig.  6.  The  marine  inversion  occurs  at 
Z  =  200  meters.  There  is  also  a  strong  layer  of  turbulence  which  occurs 
above  a  sharp  temperature  discontinuity  at  Z  -  1700  meters.   In  Fig.  8 

the  well  defined  surface  layer  extends  only  as  high  as  Z  -  30  meters. 

2 
Note  the  strong  peaks  in  CL  which  occur  at  the  dew  point  discontinuities , 

indicating  the  importance  of  water  vapor  in  atmospheric  stability.  In 

Fig.  9  we  find  very  low  levels  of  temperature  turbulence  above  the  inversion, 

2 
but  below  the  inversion  the  values  of  CT  are  very  large  considering  that 

2 
these  are  stable  conditions.  Figs.  10  and  11  show  low  values  of  CL,  near 

2 
the  surface  with  CL,  increasing  with  height  as  we  approach  maximum  tempera- 
ture at  Z  =^  400  meters.   In  this  case,  the  normal  surface  layer  equations 
are  a  very  poor  representation.   It  is  also  interesting  to  note  that  a 
nearby  profile  (Fig.  12,  taken  about  100  km  south  of  those  shown  in  Figs. 
10  and  11)  is  completely  different. 


10 


V.   CONCLUSIONS 

In  the  presence  of  a  raised  marine  inversion,  such  as  off  the  Pacific 

Coast,  a  well  mixed  turbulent  boundary  layer  is  usually  found.  The  height 

2 
dependence  of  C™  in  this  layer  will  be  well  described  by  the  standard 

surface  layer  expressions  up  to  the  inversion.   In  the  absence  of  a  strong 

raised  inversion,  such  as  is  often  the  case  off  the  Atlantic  Coast,  there 

may  be  no  well  mixed  turbulent  boundary  layer.  This  is  particularly  true 

2 
under  stable  conditions  where  the  magnitude  of  CT  found  at  the  surface 

may  in  fact  be  dominated  by  a  low  level  temperature  discontinuity.  Under 

these  conditions,  the  standard  surface  layer  equations  may  be  invalid  above 

heights  on  the  order  of  10  meters. 

Above  the  boundary  layer,  air  mass  boundaries  and  other  sources  of 

temperature,  velocity,  and  water  vapor  discontinuities  play  a  critical  role 

2 
in  the  magnitude  of  CT  .  This  is  even  more  significant  for  optical  propa- 

2 
gation  because  C,.  is  also  affected  by  water  vapor  fluctuations  (eq.  1). 

A  profile  taken  during  FLORIDA  77  (Fig.  13)  shows  that  these  layers  can 

produce  large  effects  as  high  as  5,000  meters. 
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FIGURE  CAPTIONS 

2   2 

1.  Height  dependence  of  C^  /CL,  based  on  eq.  2  for  the  U.S.  Standard 

Atmosphere. 

2 

2.  Height  dependence  of  CT  .  The  solid  line  is  the  ground  average  com- 
piled by  Hall  (1977)  for  daytime  overland  profiles.  The  dashed  line 
is  an  extrapolation  using  eq.  11  from  typical  shipboard  oceanic  sur- 
face layer  measurements. 

3.  Airborne  Research  Associates,  Inc.,  Bellanca  research  aircraft  during 
operations  off  Nova  Scotia  with  the  USNS  Hayes  in  May  of  1977. 

2 

4.  Wingtip  probe  configuration  for  measurement  of  CT  . 

2 

5.  Height  dependence  of  CL,  measured  from  the  Naval  Coastal  Systems 

Laboratory  Stage  I  off  Panama  City,  Florida  during  FLORIDA  77. 

2 

6.  Height  dependence  of  C„  for  two  profiles  during  CEWCOM  76  with  a 

simultaneous  ship  launched  radiosonde. 

2 

7.  Temperature  and  CT  profiles  off  the  New  Jersey  Coast,  22  Feb.  1977. 

The  top  of  the  haze  layer  was  1700  meters. 

2 

8.  Temper ature,dewpoint  and  CT  profile  (HAYES  77)  near  the  Nantucket 

Light  Ship  on  16  May  1977  at  1325  ADST. 

2 

9.  Temperature,  dewpoint  and  C„  profile  (HAYES  77)  near  Cape  Sable,  N.S 

2 
on  17  May  1977  at  1525  ADST.  The  x's  denote  CT  values  and  the  arrow 

denotes  the  ocean  temperature  measured  from  the  USNS  Hayes . 

2 
10.  Temperature,  dewpoint  and  C„  profile  (HAYES  77)  about  20  km  east  of 

Cape  Canso,  N.S.  on  18  May  1977  at  1625  ADST.   The  x's  denote  Cp2 

values  and  the  arrow  denotes  the  ocean  temperature  measured  from  the 

USNS  Hayes. 
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2 

11.  Temperature,  dewpoint  and  (L,  profile  (HAYES  77)  about  20  km  east 

of  Cape  Canso,  N.S.  on  18  May  1977  at  1715  ADST.  The  x»s  denote 

2 
CT  values  and  the  arrow  denotes  the  ocean  temperature  measured  from 

the  USNS  Hayes. 

2 

12.  Temperature,  dewpoint  and  Cp  profile  (HAYES  77)  about  100  km  southeast 

of  Halifax,  N.S.  on  18  May  1977  at  1750  ADST. 

2 

13.  Temperature,  dewpoint  and  CL,  profile  (FLORIDA  77)  south  of  Panama  City, 

Florida  on  19  February  1977  at  1300  CDT.  The  values  of  T  >  T  are  not 
real  but  represent  a  calibration  shift  of  the  dewpoint  instrument. 
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FIGURE   9 
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FIGURE   11 
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